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ABSTRACT 

We present a study on the effects of the intracluster medium (ICM) on the interstellar medium 
(ISM) of 10 Virgo cluster galaxies using Spitzer far- infrared (FIR) and VLA radio continuum imaging. 
Relying on the FIR-radio correlation within normal galaxies, we use our infrared data to create model 
radio maps which we compare to the observed radio images. For 6 of our sample galaxies we find 
regions along their outer edges that are highly deficient in the radio compared with our models. 
We also detect FIR emission slightly beyond the observed radio disk along these outer edges. We 
believe these observations are the signatures of ICM ram pressure. For NGC 4522 we find the radio 
deficit region to lie just exterior to a region of high radio polarization and flat radio spectral index, 
although the total 20 cm radio continuum in this region does not appear strongly enhanced. These 
characteristics seem consistent for other galaxies with radio polarization data in the literature. The 
strength of the radio deficit is inversely correlated with the time since peak pressure as inferred from 
stellar population studies and gas stripping simulations, suggesting the strength of the radio deficit 
is good indicator of the strength of the current ram pressure. We also find that galaxies having local 
radio deficits appear to have enhanced global radio fluxes. Our preferred physical picture is that the 
observed radio deficit regions arise from the ICM wind sweeping away cosmic-ray (CR) electrons and 
the associated magnetic field, thereby creating synchrotron tails as observed for some of our galaxies. 
We propose that CR particles are also re-accclcratcd by ICM-driven shocklets behind the observed 
radio deficit regions which in turn enhances the remaining radio disk brightness. The high radio 
polarization and lack of precisely coincident enhancement in the total synchrotron power for these 
regions suggests shearing, and possibly mild compression of the magnetic field, as the ICM wind drags 
and stretches the leading edge of the ISM. 

Subject headings: galaxies: clusters: individual (Virgo) — infrared: galaxies — radio continuum: 
galaxies — cosmic-rays — galaxies: interactions — galaxies: ISM 



1. INTRODUCTION 

The physical processes associated with interactions 
between the intracluster medium (ICM) and the 
interstellar medium (ISM) play critical roles driv- 
ing the evolution of spiral galaxies in clusters (e.g . 
Gunn fc Gott lll97l iLarson. Tinslev. fc Caldwell I ll98C 



Abadi. Moore fc Bower I 119991: iSchulz fc Struck I 12001 
VoUmer et al. I l2001h . Many basic effects related to 



ICM-ISM interactions (e.g., ram pressure stripping) are 
still not well understood. These effects include the 
fate of star-forming molecular clouds, the rates of trig- 
gered star formation versus gas removal, and the pos- 
sible reconfiguration of a galaxy's large-scale magnetic 
field and/or cosmic-ray (CR) particles. It is also not 
known how far from cluster centers the stripping oc- 
curs; for example , gas-deficient galax ies are found in 
cluster outskirts ( Sola nes et al. I l2001h yet it remains 
unclear if they are stripped in cluster cores or locally 
in the outskirts as seems to be the cas e in NGC 4522 
(jKennev. van Gorkom. fc Vollme71l2004h . 

If we can find clear tracers of ongoing ram pressure and 
identify galaxies in different stages of the stripping pro- 
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cess, then we can further understand the consequences of 
stripping. Since ram pressure will more easily affect lower 
density constituents of the ISM, it seems that the diffuse 
radio continuum halos of galaxies may be sensitive trac- 
ers of active ICM pressure as indicated by observations 
for a number of Virgo c luster spirals including NGC 452 2 
(|Vollmer et al. Il2004ah. NGC 4402 dCrowl et al. 1 120051 
NG C 4254 (IChvzv Ehle fc Beck 1 12007ft. and a few oth- 
ers (|Vollmer et al. I I2007D . More specifically, large re- 
gions of enhanced polarized radio continuum emission 
have been found wi t hin a number of cluster s pirals 
(|Vollmer et al. I [20071 : IChvzv. Ehle. fc Beckl [200l : the 
maxima of the polarized radio continuum within these 
galaxies are located along outer edges and thought to 
arise from external influences of the cluster environ- 
ment. However, without a good idea of the galaxy's 
unperturbed appearance in the radio, these observa- 
tions alone make it difficult to quantify the ram pres- 
sure effects. We now address this by predicting the ra- 
dio distribution expected in the absence of ram pres- 
sure by using the nearly universal correlation between 



emission of normal galaxies (e.g. 


de Jong et al. 1 


Hclou. Soifcr. & Rowan-Robinson 


1985H. 



Studies of the FIR-radio correlation with in field - 
galaxies (e.g. iHoernes. Berkhuiisen fc Xu~l 119981 : 
iHippelein et al. I 120031 : iMurphv et al. I l2006al ) have 
found that the dispersion in the FIR/radio ratios across 
a galaxy disk is as small as the dispersion of the global 
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correlation (i.e. <0.3 dex) on scales down to ^100 pc. 
While a number of other authors have studied how 
glob al FIR/radio ratios behave among cluster galaxies 
(e.g. lMiller fc Owen ||200H : iReddv fc Yun 112004ft . finding 
them to be systematically lower than those for field 
galaxies, it is only now by using Spitzer that we can 
study the influence of the cluster environment on the 
FIR/radio behavior within galaxies allowing us to better 
explain the global FIR/radio observations. 

The FIR-radio correlation is driven by the process of 
massive star formation. Massive stars both heat dust, 
yielding FIR emission, and end their lives as supernovae 
(SNe) whose remnants (SNRs), via shocks, accelerate CR 
electrons in a galaxy's magnetic field bringing about syn- 
chrotron radiation. Using a phenomenological smearing 
model, in which a galaxy's FIR map is smoothed by a 
parameterized kernel to compensate for the fact that the 
mean free path of dust heating photons is much shorter 
than t he diffusion length of CR electrons. IMurphv et al. I 
(2008) has shown for a sample of 15 non- Virgo spiral 
galaxies that the dispersion in the FIR/radio ratios on 
sub-kiloparsec scales within galaxies can be reduced by a 
factor of ~2, on average. Accordingly, it is possible to de- 
rive a good first order approximation of an undisturbed 
galaxy's non-thermal radio continuum morphology from 
its FIR image. 

Assuming that dust and gas are well mixed within 
galaxies, the FIR emission distribution provides a reason- 
able picture for the molecular gas distribution as mod- 
ulated by heating by starlight. The non-thermal radio 
continuum emission, on the other hand, is a function of 
a galaxy's magnetic field and CR electron distributions 
(i.e. the relativistic phase of the ISM). Systematic depar- 
tures from the nominal FIR/radio ratios across galaxy 
disks undergoing ICM-ISM interactions can provide in- 
sight into differential influence of these interactions on 
the gaseous and relativistic phases of the ISM. These dif- 
ferences may help to quantify the ICM pressure acting 
on the ISM of such galaxies. 

With the capabilities of current observatories the Virgo 
cluster remains an ideal laboratory allowing for many res- 
olution elements across a large number of galaxy disks, 
especially at FIR wavelengths. Coupling FIR observa- 
tions taken by the Spitzer Space Telescope with VLA 
radio continuum imaging, obtained as part of the VLA 
Imaging of Virgo in Atomic gas (VIVA; A. Chung et al. 
2008, in preparation) survey, we study how the relativis- 
tic and gaseous phases of the ISM are affected by ICM- 
ISM interactions for a sample of Virgo cluster galaxies. 
This is done using FIR Spitzer maps to predict how the 
radio morphology should appear if the galaxy were not 
experiencing an interaction; in this paper we make the 
case that significant deviations from such an appearance 
are directly related to ICM-ISM interactions. 

The paper is organized as follows: In §2 we briefly 
describe the galaxy sample and observations. We then 
present our analysis techniques for comparing the FIR 
and radio continuum imaging in §3. Results from this 
comparison are presented in §4. In §5 we compare our 
results in detail with other observational evidence of ICM 
effects for NGC 4522, a well studied galaxy experiencing 
ram-pressure stripping, and put together a physical sce- 
nario for our observations. Then, in §6, we compare and 
discuss this physical picture for the rest of the sample. 



Finally, we summarize our conclusions in §7. 

2. OBSERVATIONS AND DATA REDUCTION 
2.1. Galaxy Sample 

We present Spitzer far-infrared and VLA radio contin- 
uum imagery of 10 Virgo cluster spirals galaxies included 
in the VLA Imaging of Virgo In Atomic gas (VIVA; A. 
Chung et al. 2008, in preparation) survey. A total of 
^40 of the 53 VIVA sample galaxies are included in the 
Spitzer Survey of Virgo (SPITSOV; J. D. P. Kenney et al. 
2008, in preparation) imaging program. The sub-sample 
of 10 galaxies were the first few observed with Spitzer 
that also had high quality radio maps as well as FIR 
and radio emission detected out about as far as the H 1 
extent. Because we are interested in studying the prop- 
erties of the radio continuum halo to determine if it is 
interacting with the ICM, we use this criterion to ensure 
we are detecting infrared and radio continuum emission 
near the edge of the ISM. 

Basic data for each galaxy can be found in Table [TJ 
Galaxy diameters (-D25), morphological types, position 
angles (PA), and absolute B-band magnitudes, assum- 
ing a distance of 16.6 Mpc to the Vi rgo cluster center 
(jShaplev. Fabbiano. fc Eskridge Il200lh . were taken from 
the Third Reference Catalo g of Bright Galaxies (RC3; 
Ide Vaucouleurs. et al. Ill99"lt ). 

We also include distanc es to the Virgo clust er center 
(12 h 30 m 47 s , +12°20'13"; lEbeling et all fl998h . in de- 
grees, for which the sample ranges between 1.3 — 7.2° 
(i.e. spanning a factor of >5). Since the Hubble mor- 
phological classification appears to have limited meaning 
for cluster spirals (e.g. stripped Sc's te nd to be classi- 
fied as Sa's; lKoopmann fc Kennev 111998ft. we include th e 
Ha star formation type of lKoopmann fc Kennev I (2004) . 
Type-independent H I deficiencies, taken from A. Chung 
et al. (2008, in preparation), are included in Table Q] and 
have a typical uncertainty of 0.2 dex arising from the 
uncertainty in morphological classifications. The present 
sample spans a factor of ~50 in H 1 deficiency and nearly 
the entire range in Ha star formation type. 

To determine if the FIR and radio properties of these 
cluster galaxies differ systematically from those in the 
field, we compare with 6 non- Virgo galaxies (NGC 2403, 
NGC 3031, NGC 3627, NGC 4631, NGC 5194, and 
NGC 6946), chosen to act as a control sample. These 
galaxies were selected to be less distant than Virgo, span 
a range in spiral type and inclination, and are a sub- 
sa mple of the galaxies in cluded the FIR and radio study 
of IMurphv et al. 1 (|2008f ); this investigation describes the 
physical link between the FIR and radio properties of 
normal galaxy disks on kiloparsec scales. Information 
on the observations and general characteristics of these 
galaxies can be found there. 

2.2. Spitzer MIPS Imaging 

Spitzer imaging was carried out for each galaxy using 
the Multiband Ima ging Photometer for Spitzer (MIPS; 
iRieke. et all 120041) . The strategy of the imaging cam- 
paign was based on that o f the Spitzer Infrared N earby 
Galaxies Survey (SINGS; iKennicutt et al. I 12003ft with 
the only difference being a factor of 2 increase in expo- 
sure times to better detect diffuse emission arising in the 
outer regions of the Virgo sample galaxies. The MIPS 
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TABLE 1 
Basic Galaxy Data 





R.A. 


Decl. 








M B 


V r 


D 25 




PA 


87 


H I Def. 


Galaxy 


(J2000) 


(J2000) 


Type 


Ha Type 


Nuc. 


mag 


(km s~ v ) 


(arcmin) 


b/a 


(°) 


(°) 


(dex) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


(13) 


NGC 4254 


12 18 49.4 


+14 25 07 


SAc 


N 




-20.66 


2407 


5.4 


0.87 





3.6 


-0.10 


NGC 4321 


12 22 55.2 


+15 49 23 


SABbc 


T/N 
T/N a 


L/H 


-21.05 


1571 


7.4 


0.85 


30 


4.0 


0.35 


NGC 4330 


12 23 16.5 


+11 22 06 


Scd 


-18.01 


1565 


4.5 


0.20 


59 


2.1 


0.80 


NGC 4388 


12 25 47.0 


+12 39 42 


SAb 


T/N [s] a 


Sy2 


-19.34 


2524 


5.6 


0.23 


92 


1.3 


1.16 


NGC 4396 


12 25 59.3 


+15 40 19 


SAd 


N a 




-18.04 


-128 


3.3 


0.29 


125 


3.5 


0.30 


NGC 4402 


12 26 07.9 


+13 06 46 


Sb 


T/N a 




-18.55 


232 


3.9 


0.28 


90 


1.4 


0.74 


NGC 4522 


12 33 40.0 


+09 10 30 


SBcd 


T/N [s] 




-18.11 


2307 


3.7 


0.27 


33 


3.3 


0.86 


NGC 4569 


12 36 50.1 


+13 09 48 


SABab 


T/N [s] 


L/Sy 


-20.84 


-235 


9.6 


0.46 


2.3 


1.7 


1.47 


NGC 4579 


12 37 44.2 


+11 49 11 


SABb 


T/N 


L/Syl.9 


-20.62 


1519 


5.9 


0.79 


95 


1.8 


0.95 


NGC 4580 


12 37 48.4 


+05 22 09 


SABa pec 


T/N [s] 




-19.27 


1034 


2.1 


0.78 


165 


7.2 


1.53 



Note. — Col. (1): ID. Col. (2): The ri ght ascension in the J2000.0 e poch. Col. (3): The declination in the J2000.0 epoch. Col. (4): 
RC3 type. Col. (5): Star formation class of Koopmann & Kcrmcy (2004): N: Normal; A: Anemic; E: Enhanced; T/A: Truncated/Anemic; 
T/C: Truncated/Compact; T/E: Truncated/Enhanced; T/N: Truncated/Normal; T/N [s] Truncated/Normal (Severe); a : classified for this 
paper. Col. (6): Nuclear type: H: H II region; L: LINER; Sy: Seyfert (1, 1.9, 2); SB: Starburst. Col. (7): Absolute S-band magnitudes. 
Col. (8): Heliocentric velocity. Col. (9): Major axis diameters. Col. (10): Semi- minor/semi- major axis ratio. Col. (11): Position 
Angle in degrees. Col. (12): Distance from Virgo cluster center (12 h 30 m 47 s . +12°20'13"; lEbeling et al. Ill99ll) in degrees. Col. (13): 
Type-independent H I deficiencies, having a typical uncertainty of 0.2 dex arising from the uncertainty in morphological classifications, 
taken from A. Chung et al. (2008, in preparation) (deficiencies logarithmically increase with number). 



data were reduced as in iMazzarella et al.~l (|2008f ) with 
the following modifications. Since most of the galax- 
ies were observed twice, the data had to be combined 
into a single mosaic of the two observations. For the 
70 and 160 /im data this was a straightforward pro- 
cess, combining the filtered basic calibrated data frames 
(BCDs), which have near zero residual background, from 
the two observations using the Mosaicking and Point 
Source Extraction (MOPEX) software package. The tar- 
get gabxy_was_jri^kedduring the 70 /zm filtering, as 
in IMazzarella et al.~l ([2008) , while the pipeline filtered 
BCDs (FBCDs) were used for 160 (im observations. The 
individual 24 /im observations were first reduced sepa- 
rately in order to measure the background level. The 
background level was then subtracted from each BCD, 
and a combined mosaic was constructed using MOPEX 
as above. Although the two observations for a given 
galaxy were often performed in close succession, this was 
not always the case. Measuring and subtracting the two 
backgrounds was the most reliable method of placing the 
two observations on a common scale for mosaicking. 

For galaxies which were observed once at the time of 
this writing, preliminary results were obtained by com- 
bining the BCDs (as above, filtered in the case of 70 /im 
and 160 /im) using MOPEX. The resulting 160 /im mo- 
saics have small regularly spaced coverage gaps at the 
locations of bad pixels on the array. This is a conse- 
quence of the single redundancy of 160 /im scan mode 
observations, although such occurrences did not affect 
the field of view of our galaxies in any observation. The 
final calibration uncertainties of 24, 70, and 160 /im data 
are 5, 10, and 15% respectively. 

For galaxies in common with the VIVA survey and 
SINGS, we use MIPS data included in the SINGS 
data release 5 (DR5) (see Tabled]). These data were 
processed using the MIPS Data Analysis Tool (DAT; 
iGordon et ani2005T ) version 3.06. Additional steps be- 
yond the standa rd reduction pro c edure o f the MIPS DAT 
are described in iMurphv et al. 1 (|2006a|) . The full width 
at half maximum (FWHM) of the MIPS 24, 70, and 
160 /im point spread functions (PSFs) are ~5'.'7, 17" 



and 38" while the final calibration uncertainties are ~2, 
5, and 9%, respectively. The mean RMS noise values 
of the 24, 70, and 160 /im data for the SINGS galaxies 
are ~0.048, 0.45, and 0.57 MJy sr^ 1 while our deeper 
imaging campaign reaches mean RMS noise values of 
~V2 better. In some cases (e.g. NGC 4254, NGC 4321, 
NGC 4579), bright nuclei created negative latent images 
in the 70 /im data resulting in bright and dark streaking 
on opposite sides of the nuclei after background subtrac- 
tion. While visible in some of the images, these streaks 
are generally low amplitude and should not significantly 
affect the analysis. 

2.3. Radio Imaging 

The 1.4 GHz radio continuum maps were created from 
the line-free channels of H I data cubes collected as part 
of of the VIVA survey; a detailed description of the VLA 
observations along with the H I reductions and associated 
data products can be found in A. Chung et al. (2008, in 
preparation). In order to reduce the effects of interfer- 
ing sources, which causes substantial sidelobes especially 
at low frequencies, we have used the AIPS procedure 
PEELR. It iteratively attempts to calibrate on multiple 
fields around bright continuum sources (self-calibration) , 
subtract those fields from the self-calibrated data, undo 
the field-specific calibration from the residual data, and 
it finally restores all fields to the residual data. The 
quality of our continuum data varies depending on the 
number of line free channels of individual galaxies. The 
same weighting scheme (R0BUST= 1) used for the the H I 
maps was also used here for the continuum maps. The 
RMS of each 20 cm radio continuum map is given in 
Table H The mean RMS noise value of the VIVA radio 
maps is ~0.008 MJy sr _1 (0.1 mJy bnr 1 ) and the range 
in RMS noise values spans a factor of ^3. The FWHM 
of the CLEAN beam major and minor axes, -B ma j and 
B m in respectively, are also given for each galaxy in Table 
[2 the major axes range from ~17" — 45". 

In the case of NGC 4569, a SINGS galaxy observed 
as part of the Westerb ork Synthesis Rad io Telescope 
(WSRT)-SINGS survey (|Braun et al. 112001 . we use the 
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TABLE 2 

Observation Summary 





MIPS 




20 cm RMS 


70 /im RMS 


Galaxy 


Summary (arcsecond) 


MJy sr _1 


MJy sr _1 


(1) 


(2) 


(3) 


(4) 


(5) 




Virgo (Cluster) Sample 


NGC 4254 


s 


34.9 x 23.3 


0.0066 


0.45 


NGC 4321 


s 


27.0 x 22.3 


0.0066 


0.24 


NGC 4330 


1 


19.4 x 15.4 


0.0045 


0.30 


NGC 4388 


2 


17.3 x 15.1 


0.014 


0.37 


NGC 4396 


2 


20.5 x 20.2 


0.0047 


0.18 


NGC 4402 


1 


17.1 x 15.2 


0.0097 


0.54 


NGC 4522 


2 


19.0 x 15.3 


0.0066 


0.26 


NGC 4569 6 


s 


45.0 x 45.0 


0.0015 


0.22 


NGC 4579 


s 


43.2 x 35.2 


0.0048 


0.26 


NGC 4580 


1 


17.5 x 16.6 


0.0094 


0.50 


Field (Control) Sample c 


NGC 2403 


s 


92.7 x 92.7 


0.0020 


0.13 


NGC 3031 


s 


82.2 x 82.2 


0.0028 


0.13 


NGC 3627 


s 


31.8 x 31.8 


0.0039 


0.45 


NGC 4631 


s 


39.0 x 39.0 


0.0023 


0.24 


NGC 5194 


s 


38.4 x 38.4 


0.0015 


0.24 


NGC 6946 


s 


43.9 x 43.9 


0.0036 


0.35 



Note. — Col. (1): ID. Col. (2): Spitzer MIPS Data: Si Observed 
as part of SINGS; 1: Only a single AOR included; 2: Both AORs in- 
cluded. Col. (3): Major and minor axes of the original radio CLEAN 
beam used in the VIVA radio continuum maps (A. Chung et al. 2008, 
in preparation). Col. (4): 20 cm RMS. Col. (5): 70 (im RMS. 
a : The FWHM of the PSF at 70 /im is ~17" for each galaxy; B maj 
therefore sets the resolution for the analysis of each galaxy disk. 
b : For NGC 4 569 w e use the WSRT-SINGS radio continuum map of 
IBraun et al. I 1120071 ) rather than the VIVA radio continuum data be- 
cause the WSRT-SINGS map is signifi cantly more sen sitive. 
c : These data were taken from Murphy et al. (2008;); the resolution 
indicates that for which the galaxy's maps were smoothed to simulate 
their appearance at the distance of Virgo. It is at this resolution for 
which the RMS noise was measured. 



WSRT radio map in our analysis since it is an order of 
magnitude more sensitive than what was achieved by the 
VIVA survey for this object. 

2.4. Image Registration and Resolution Matching 

To perform an accurate comparative analysis between 
the MIPS and radio data we match the resolution of the 
final calibrated images using the MIPS PSF. After crop- 
ping each set of galaxy images to a common field of view 
we CLEANed the radio data and convolved the resulting 
CLEAN components with a model of the MIPS 70 /im 
PSF. We restore with the MIPS beam first because it has 
a non-Gaussian shape with significant side-lobes. 

Since the FWHM of the original radio CLEAN beam 
-Bmaj was in all cases larger than the FWHM of the MIPS 
70 /im PSF, we convolved the 70 /im and 20 cm data with 
an additional Gaussian beam to account for this differ- 
ence. These images were then re-gridded to a common 
pixel scale and properly aligned. The observed 70 /im 
and resolution-matched 20 cm maps for each galaxy are 
presented in the first and second columns of Figure [TJ 
The effective resolution used in the comparison of the 
infrared and radio data is given in Table O 

2.5. Global Flux Densities and FIR/Radio Ratios 

In Table we give the global MIPS and 20 cm flux 
densities. We include extended source aperture correc- 



tions for the MIPS fluxes b y ado pting the median val- 
ues derived bv iDale et aT~l (j2007f ) which are 1.01, 1.04, 
and 1.10 at 24, 70, and 160 /im respectively. Because 
our sample does not include dwarf galaxies, for which 
the corrections were found to be large, simply taking the 
median value should be adequate. Since we are using 
the WSRT-SINGS radio continuum map of NGC 4569, 
"we also w ork with the g l obal fl ux density at 22 cm rc- 
- ported by IBraun et al. I (|2007ft which has been scaled 
to what is expected at 20 cm assuming a synchrotron 
power law of the form S v oc v~ a with a radio spec- 
tral index a = 0.8 (i.e. typical for s tar-forming galaxies; 
iNiklas. Klein, fc WielebinskVI[l997h . This was also done 
for the 22 cm flux densities of the WSRT-SINGS field 
galaxies where appropriate. 

We also provide measures of the total infrared 
flux (TIR; 3 - 1100 /im) along with logarith- 
"mic far-infrared (FIR; 42 — 122 /im)-to-radio ratios 
- in the form of the commonly used parameter q 
(Hcl ou. Soifer. fc Rowan- Robinson 111985ft such that, 



q = log 



F 



FIR 



3.75 x 10 12 Wm- 2 



log 



/ Si 



4 GHz 



VWnr 2 Hz- 



The TIR fluxes were measured using a linear combi- 
nation of the 3 MIP S bands given by Equation 4 of 
IDale fc Heloul (|2002fl . The corresponding FIR fractions 
were then calculated usin g the same spec tral energy dis- 
tribution (SED) models oflDale fc Heloul (12002ft as given 
by Equation 3 of lMurphv et al. I (|2008ft . i.e.. 



F 



FIR 



TIR 



3 

£ 

i=0 



6 log 



UJO A™) 
/„(160 Mm) 



(2) 



where [£0,6,6,6] = [0.5158,0.1619,-0.3158,-0.1418] 
and -0.65 < log /„(70 /im)//„(160 /mi) < 0.54. 

3. IDENTIFYING ICM RELATED EFFECTS IN 
THE RADIO MORPHOLOGIES 

3.1. Radio and Infrared Morphologies 

A comparison of the 70 /im and 20 cm images in the 
first two columns of Figure [T] shows that the radio contin- 
uum images for some of the galaxies are quite disturbed, 
appearing asymmetric compared to the infrared images 
(e.g. NGC 4330 and NGC 4522). In some cases we find 
regions along edges having FIR emission without corre- 
sponding radio emission (these are discussed further in 
§6.4); regions opposite of these edges often have radio 
continuum tails without FIR counterparts. To illustrate 
this we plot pixel identification maps for NGC 4402 and 
NGC 4522 in Figure^ regions of cyan indicate that both 
the radio and 70 /im emission was detected at the 3-c 
RMS level while region of red and yellow indicate areas 
where only the radio or only the 70 /im emission has been 
detected at the 3-c RMS level. 

The FIR surface brightness depends on the ISM 
density and heating intensity by starlight. Assuming 
that the FIR dust emission is, to first order, tracing 
the distribution of molecular gas within galaxies (e.g. 
Murg iaet al. 112005ft . the observed differences in the ra- 
dio continuum and FIR emission distributions suggest 
that the relativistic ISM is perturbed in cluster galaxies 
more easily than molecular material. We would like now 
to interpret and quantify these effects. 
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TABLE 3 

Global Flux Densities and Derived Parameters 





S„(20 cm) 


/„(24 Mm) 


/„(70 fim) 


/„(160 Atm) 






Galaxy 


(Jy) 


(Jy) 


(Jy) 


(Jy) 


(10- 13 W m" 2 ) 










Virj 


;o (Cluster) Sample 






NGC 4254 


0.439 


4.17 


45.0 


147 


60.0 


2.13 


NGC 4321 


0.273 


3.44 


36.8 


143 


54.9 


2.26 


NGC 4330 


0.020 


0.107 


1.05 


5.65 


1.97 


1.92 


NGC 4388 


0.164 


2.64 


10.6 


20.0 


13.7 


2.00 


NGC 4396 


0.022 


0.132 


1.74 


5.79 


2.29 


2.00 


NGC 4402 


0.067 


0.689 


7.52 


25.5 


10.2 


2.17 


NGC 4522 


0.027 


0.217 


1.43 


6.60 


2.55 


1.91 


NGC 4569 


0.157 a 


1.41 


11.2 


41.5 


16.8 


2.00 


NGC 4579 


0.167 


0.816 


9.34 


42.8 


15.4 


1.89 


NGC 4580 


0.005 


0.148 


1.83 


7.75 


2.84 


2.73 


Field (Control) Sample 6 


NGC 2403 


0.33 a 


5.84 


86.4 


246 


102 


2.50 


NGC 3031 


0.55 a 


5.09 


85.2 


360 


129 


2.31 


NGC 3627 


0.46 a 


7.42 


92.6 


230 


103 


2.38 


NGC 4631 


1.20 a 


8.15 


130 


290 


132 


2.09 


NGC 5194 


1.15 a ' c 


12.7 


147 


495 


198 


2.22 


NGC 6946 


1.57° 


20.4 


207 


503 


234 


2.21 



Note. — a : Scaled from the 22 cm flux density reported bv lBraun et al. I lf2007l) assuming 
a radio spectral index a = 0.8 (see jJSHJ. 



b : D ata taken from Tab le 2 of lMurphv et al.1 fl2008T ). where MIPS flux densities were taken 
from lDale et al.H 1120071 ). q values have been adjusted since here we are using 20 cm radio 
flux densities rather than 22 cm flux densities. 

c : Corrected for flux contributions from companion galaxy NGC 5195. 

factor of ^2 using an image-smearing technique which 
approximates a galaxy radio image as a smoother ver- 
sion of its FIR image due to the diffusion of CR elec- 
trons. While the star-forming disks of cluster galaxies 
are often truncated due to gas stripping events, the im- 
ages of the remaining 70 /im disks for the entire sample do 
not appear to be otherwise more strongly disturbed than 
galaxies in the field; we do not observe any obvious sharp 
edges, tails, or evidence for strong extraplanar emission 
at 70 /x m. Therefore, we apply the image-smearing tech- 
nique of lMurphv et al. I (|2008f ) (using a single smoothing 
function) to create models for the expected radio distri- 
bution of a galaxy assuming that the FIR/radio ratio 
map is like that of a field galaxy; the modeled radio con- 
tinuum images are denoted as i? mo d- 

Since the image-smearing description is applicable 
purely for non-thermal radio emission, we include an es- 
timate of the thermal radio emission i? x derived from 
the 24 jim f lux density distribution (see iMurphv et al. I 
I2006U 120081 for details) such that, 

m ~ 7.93x10-3 f^_r 5 (^r - 1 , 

(3) 

where we have assumed an average H n region tempera- 
ture of T = 10 4 K. 

We compare the observed radio continuum images with 
the model radio continuum images, i? m od, such that 

Rmod = R-T + Qnt(«(0 * J )> ( 4 ) 

where i?x is an estimate of the thermal radio emis- 
sion distribution using Equation [5J J is the 70 /im im- 
age, * denotes a convolution with an exponential ker- 
nel k having scale-length I, and Qnt is the global 
/„(70 /im) / S u (20 cm) flux density ratio after correcting 



TABLE 4 
Image-Smearing 
Parameters 



I <S> 
Galaxy (kpc) (dex) 



Virgo (Cluster) Sample 



NGC 4254 


2.3 


0.927 


NGC 4321 


0.7 


0.253 


NGC 4330 


0.7 


0.381 


NGC 4388 


0.2 


0.078 


NGC 4396 


0.9 


0.668 


NGC 4402 


0.6 


0.581 


NGC 4522 


0.8 


0.835 


NGC 4569 


0.8 


0.933 


NGC 4579 


0.0 


0.000 


NGC 4580 


0.8 


0.762 


Field (Control) Sample" 


NGC 2403 


1.0 


0.664 


NGC 3031 


1.8 


0.464 


NGC 3627 


0.5 


0.541 


NGC 4631 


0.6 


0.768 


NGC 5194 


0.5 


0.343 


NGC 6946 


0.6 


0.590 



N PTE. — a : Values tak en 
from IMurphv et al. I (2008). 



3.2. Modeled Radio Continuum Maps 

Detailed studies of the FIR-radio correlation within 
nearby field galaxies (e.g. IMurphv et aFI l2006al I2008D 
have shown that the dispersion in the FIR/radio ratios 
on >0.1 — 1 kpc scales within galaxy disks is similar to 
the dispersion in the global FIR-radio correlation (i.e. 
<0.3 dex). This strong correlation between the FIR and 
radio images of spirals has been found to improve by a 
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Fig. 1. — Each galaxy's 70 jim and 22 cm images are displayed in columns 1 and 2, respectively. Streaking artifacts in the 70 fim data 
of NGC 4254, NGC 4321, and NGC 4579 are easily identifiable. The stretch of gray-scale and contours ranges, logarithmically, from the 
3-cr RMS level of the bac kgro und to the maximum surface brightness of the galaxy. In column 3 the logarithm of the ratio of observed to 
predicted radio flux (see i|3,2H is shown with a stretch ranging from 10% (dark) to a factor of 3 (light). The radio deficit and excess regions 
(see £|3.3| I are given in column 4 along with a vector showing the direction to cluster center. Overlaid on each deficit (black) and excess 
(cyan) region map are H I contours; the H I contours begin at the ~3-cr RMS level and increase logarithmically. Using a dashed line we 
overlay the extent of the optical disk on the 70 (im, 22 cm, and H I maps. A 5 kpc scale bar is also drawn in each of the 70 (im image 
panels, assuming a distance of 16.6 Mpc to the Virgo cluster center. 



the global 20 cm flux density for thermal radio emis- residuals betwee n the radio and smoothed in frared maps 
sion using the galaxy's 24 fim flux density and Equation as described by (|Murphv et al. I feOOGa. 2008f) ; only areas 
[3J The scale-length I is determined by minimizing the of the galaxy disk which are detected above the 3-cr RMS 
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Fig. 1. — Continued 



noise in both of the observed radio and infrared images 
are considered. Thus, the residuals are calculated over 
the non-disturbed regions of the radio disks. 

To quantify the improvement in the spatial correla- 
tion between the radio and smoothed 70 /xm images we 
measure $ which is the logarithmic decrement from the 
squared sum of the residuals between the observed (un- 
smoothed) radio and infrared images to the minimum 



residuals which defines the best-fit scale-length. The 
values of the best-fit scale-lengths used for creating the 
modeled radio maps along with the improvement in the 
residuals between the two maps, as defined by <E>, are 
given in Table [4j 

We find that the the image-smearing improves the dis- 
persion between the infrared and radio maps by a factor 
of ~2, on average; this is the same quantitative improve- 
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Fig. 2. — Pixel identification maps indicating regions of the 
galaxy for which detections were made at the 3-cr RMS level. Cyan 
regions indicate areas where both the radio and 70 /an emission 
were detected at the 3-cr RMS level while red and yellow regions 
identify where only the radio or 70 jim emission has been detected 
at the 3-CT RMS level. 

ment quoted bv lMurphv et al. 1 (|2008f) using largely field 
galaxies. We also find best-fit scale-lengths of 0.2 kpc for 
NGC 4388 and 0.0 kpc for NGC 4579, which likely result 
from their bright Seyfert-type nuclei. 

3.3. Deviations from Expectations for Field Galaxies 

Our aim is to determine if differences between the ob- 
served and modeled radio continuum images arise from 
ICM-ISM interactions and, if so, whether quantifying 
and comparing such differences within our sample can 
help improve our current physical picture of such inter- 
actions. We therefore create ratio maps between the ob- 
served and modeled radio maps in the following manner. 
We divide the observed radio map by the modeled ra- 
dio map, keeping only those pixels (r^) detected at >5-<7 
RMS level of either map. We also construct a map giving 
the associated uncertainty at each pixel of the ratio map 
(c n .RMs) by numerically propagating the uncertainties 
of the input images as measured by the 1-a RMS of their 



maps. Any background radio galaxies or artifacts aris- 
ing from streaking in the 70 /zm data (e.g. NGC 4254, 
NGC 4321, and NGC 4579) outside the galaxy disk are 
removed. The ratio maps are displayed in column 3 of 
Figure [TJ 

Next, we identify regions for which the ratio of ob- 
served to modeled radio flux density departs significantly 
from unity. To do this we compute the interquartile stan- 
dard deviation of the ratio map which is the standard de- 
viation computed after removing values included in the 
first and fourth quartiles; we denote this value as oy^isk 
which, for a Gaussian distribution, is equal to ~^ the 
population standard deviation. Radio deficit regions are 
constructed using pixels which simultaneously obey the 
following two conditions: {|1 — fil/oy^RMS > 2.5} and 
{(1 — 7"j)/oy iC iisk > 2.5}. The first condition selects pixels 
which are significantly deviant with respect to the RMS 
noise of the maps while the second condition selects pix- 
els which are significantly deviant with respect to the in- 
trinsic dispersion of the ratio map within the galaxy disk. 
Similarly, radio excess regions are constructed from all 
those pixels in the ratio map which simultaneously obey 
the conditions such that: {|1 — ri\/a ri! RMS > 2.5} and 
{(n - l)/oy !disk > 2.5}. 

Radio excesses appear largely coincident with regions 
associated with synchrotron tails like those observed in 
NGC 4330, NGC 4402 and NGC 4522. In the cases of 
NGC 4524, NGC 4321, and NGC 4579, the radio excess 
regions appear to be possibly associated with the streak- 
ing in the 70 /jm data. On the other hand, for most of 
our sample galaxies we detect radio deficit regions which 
are located along a single edge of their disks. The radio 
deficit regions are generally found opposite radio excess 
regions that are associated with synchrotron tails; this 
is illustrated in Figure [3] for NGC 4330, NGC 4402, and 
NGC 4522. NGC 4254 is an exception, for which an 
excess region is also found interior to the deficit region. 
Since these synchrotron tails are likely identifying the di- 
rection of the ICM wind, we focus our attention toward 
the radio deficit regions as they are probing directly the 
most intense effects of the ongoing ICM-ISM interactions. 

We define the maximum ratio, or fractional deficit, 
used for identifying radio deficient pixels as r^°[ = 1 — 
2-5cr I% disk and find < rjjj^ >~ 0.5 ±0.1 among the sample 
galaxies. We note that the lack of high-quality radio con- 
tinuum data for the entire sample is a limiting factor for 
this study. The large range in the RMS noise of the ra- 
dio maps (see Table [2]) prohibits a proper comparison be- 
tween the deficit regions for each galaxy. However, since 
there is such small scatter in r~fj., we choose to use a con- 
stant value of rj?ut = 0.5, or 50% of the expected flux, for 
creating the deficit region maps of all galaxies; this allows 
for the most reliable inter-galaxy comparison. Doing the 
same for excess regions we find < r^jr >^ 1.3 ± 0.1 and 
choose a constant value of r^l = 1.3, or 1.3 times the 
expected flux, for creating the excess region maps of all 
galaxies. The radio continuum deficit (black) and excess 
(cyan) regions are given in column 4 of Figure [1] and are 
overlaid by the contours of the VIVA H I images. 

To compare quantitatively the radio deficiencies among 
the sample galaxies we define two parameters. First, we 
measure the ratio of the deficit region area (Adei) to that 
of the total galaxy area given by the size of our ratio maps 
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NGC 4330 NGC 4402 NGC 4522 

Fig. 3. — The radio deficit regions of NGC 4330, NGC 4402, and NGC 4522 with radio continuum contours. The radio continuum 
contours begin at the 3-cr RMS level and increase logarithmically. 



(^4disk) which are determined by the regions of either 
significant FIR or radio emission. Next, we calculate the 
severity of the deficit as measured by the parameter 



oohs\ 
&v Jdei 



which measures the difference between the observed and 
modeled radio flux densities within the radio deficit re- 
gion normalized by the global radio flux density of the 
galaxy (i.e. the flux densities given in Table [3]). These 
values along with the flux densities and areas of the 
deficit regions are given in Table O 

4. RESULTS COMPARING OBSERVED AND 
PREDICTED RADIO DISTRIBUTIONS 

4.1. The Ratio Maps 

The ratio maps in column 3 of Figure [1] exhibit gra- 
dients near one edge of the galaxy disks for most of the 
sample; the lowest ratios lie preferentially along these 
edges. A comparison of these gradients with our beam 
demonstrates that they are in fact resolved. Further- 
more, these gradients are opposite to what is expected 
if we were not smoothing the infrared images enough 
to estimate properly the CR-electron diffusion lengths 
since, as observed, the radio disks of galaxies are gen- 
erally found to be characterize d by a larger scale-length 
than their infrared disks (e.g. iBicav. Helou. fe Condon I 
[LOS! iBicav fc Helou 111990 : iMarsh fc Helou 111995ft . Un- 
like the rest of the sample, we do not identify statis- 
tically significant radio deficiencies within the disks of 
NGC 4321 and NGC 4580. 

We do detect radio deficit regions in NGC 4396 
and NGC 4579, however these regions may arise from 
non-ICM related effects. The radio deficit region for 
NGC 4396 is detected to the north next to a number 
of nearby background galaxies. It is unclear if the ob- 
served radio deficit region partly due to a confusion ar- 
tifact resulting from the background galaxies, or if it is 
real. Unlike the infrared map, the radio continuum and 
H I images appear to have a sharp edge coincident with 
the radio deficit region; however, this region is not op- 
posite the location of extended H I emission. In the case 
of NGC 4579 we find a radio deficit region interior to 
the edge of the H I disk. The region is consistent with a 
streaking artifact in the 70 /xm image resulting from the 
galaxy's bright, Seyfert-type nucleus. 

For the remaining six galaxies (i.e. 60% of the sample), 
regions identified to be radio deficient are always at an 



edge of the disk and opposite of any identified H I tails. 
For example, H I tails are clearly identified by the con- 
tour overlay on the radio continuum deficit region maps 
of NGC 4330 and NGC 4522 in column 4 of Figure ED 
H I tails are also observed for, and opposite to the ra- 
dio deficit regions of NGC 4254, NGC 4388, NGC 4402, 
and NGC 4569, however they are not easily visible in 
Figure [1] either because the VIVA data is not sensitive 
enough to pick up the low level emission (i.e. NGC 45254 
and NGC 4388) or because of smoothing the data for 
resolution matching (i.e. NGC 4402 and NGC 4569). 
The radio continuum deficit regions also appear oppo- 
site any observed radio continuum tails which is clearly 
illustrated for NGC 4330, NGC 4402, and NGC 4522 in 
Figure [3] The combination of these observations sug- 
gests that the radio deficit regions likely arise from the 
same gravitational or gasdynamical effects which have 
displaced the galaxy's H I gas and relativistic plasma. 
Since we think that the H I tails are created by the ICM 
wind, we therefore believe that the radio deficit regions 
identify the zone in which the ICM wind is actively work- 
ing on each galaxy's ISM. 

For those sample galaxies with radio deficit regions, the 
fractional area of the disk defined to be radio deficient 
ranges from ^5 to 30%. We also find that the difference 
between the observed and modeled flux densities within 
the deficit region, normalized by the global flux density 
(i.e the parameter T), range from ~1 to 15%. Not sur- 
prisingly, we find that deficit area ratio increases with 
increasing T (see Table [5]) . 

4.2. Comparison with Field Galaxies 

To determine whether the identification of the radio 
continuum deficit regions is unique to these cluster galax- 
ies we repeat the above analysis on six non- Virgo galax- 
ies observed as part of SINGS: NGC 2403, NGC 3031, 
NGC 3627, NGC 4631, NGC 5194, and NGC 6946. A 
detailed comparison of the FIR and radio properties of 
these galaxies, along wi th the rest of the WSR T-SINGS 
sam ple, can be found in IMurphv et al. I (|2008D . We note 
that IMurphv et al. I (|2008t ) do not report finding asym- 
metric gradients in the FIR/radio ratio maps for their 
sample of mostly field galaxies, suggesting that the FIR 
and radio properties of these Virgo galaxies do differ sys- 
tematically from those in the field. To test this further, 
we repeat the exercise of looking for radio deficit regions 
among these six galaxies. 

The smoothing scale-lengths used to create the mod- 
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eled radio maps are given in Table |4j However, before 
constructing model radio and ratio maps, the infrared 
and radio images of these galaxies were smoothed appro- 
priately to simulate their appearance at the distance of 
Virgo. This ensures distance effects, such as differences 
in the physical resolution of each disk, will not introduce 
any biasses into the comparison. The final beam sizes, 
along with the corresponding noise in the radio and in- 
frared maps, are given in Table [H The ratio maps of 
each galaxy are given in Figure 2] using the same color 
scale as for the ratio maps of the Virgo galaxies in the 
third column of Figure [1] 

Unlike some of the ratio maps for the Virgo sample 
(e.g. NGC 4402 and NGC 4522) we do not find strong 
gradients near the edges of the disks for the non- Virgo 
galaxies. We also note that it was not necessary to 
show maps of radio continuum deficit regions for the 
WSRT-SINGS galaxies since, except for two cases, ra- 
dio deficits were not detected. The first exception is 
NGC 3031 where we find a few pixels (i.e. <0.8% of 
the total area) to be radio deficient just surrounding its 
nucleus; this likely arises from its bright AGN. The other 
exception is NGC 3627, in which the north-west region of 
the galaxy is significantly radio deficient (i.e. T = 0.011 
and Ajof/Ajisk = 0.08 as noted in Table [5]); among the 
Virgo galaxies with radio deficit regions detected, only 
NGC 4524 has a lower T value. The radi o deficit region 
is loc ated opposite of its observed H I tail (|Havnes et al. I 
Il979t ) and just north of a n area observed to have highly 
polarized radio emission (jSoida et al. iflOOOL [200lh . The 
consistency in the location of the deficit region with other 
tracers for ongoing ICM-ISM interactions, combined with 
the fact that NGC 3627 is a member of the Leo Triplet, 
suggests that the radio deficit region may arise from an 
interaction within the group. The control sample behav- 
ior clearly demonstrates that the radio deficiency mor- 
phology observed in Virgo galaxies is singular and most 
naturally explained as a result of ICM-ISM interactions. 

5. NGC 4522: A WELL STUDIED CASE OF RAM 
PRESSURE STRIPPING 

5.1. Atomic and Ionized Gas 

The Virgo spiral NGC 4522 is one of the clear- 
est cases where the effects of ram pressure, result- 
ing from the galaxy's rap id motion through the ICM , 
are directly observed (e.g. iKenney fe Koopmann 1 119991: 
Kennev. van Gorkom. fe Vo llmer ' I2004t lVoll~ et al. I 
2004al ICrowl fc Kennev I I2006D . A combination of ob- 
servational data in cluding asymmetric and ext raplanar 
emission in the Ha ( Kennev & Koopmannl ll999l ) and H I 
([Kennev. van Gorkom. fc Vollmer Il2004h lines are highly 
suggestive of such ICM-ISM interactions. While the ICM 
wind has stripped away atomic gas, apparent by the ex- 
traplanar H I gas, the extraplanar Ha emission is thought 
to arise from extraplanar H II regions; so it appears that 
virtually the entire ISM has been stripped away from the 
outer disk, including molecular clouds. The stripping ap- 
pears to have occurred very recently. Both the stellar 
population st udy of the stellar disk be yond the gas trun- 
cation radius (|Crowl fc Kennev1l2006f ) and a comparison 
of the ob served gas distribution and kinematics with sim- 
ulations (jVollmer et al. I [20061 ) indicate that the galaxy 
was stripped within the last 50 — 100 Myr. Thus, it is 



likely that the current ram pressure is still significant. 

The H I contour overlay on the radio deficit region 
of NGC 4522 in Figure [1] shows that the deficit region 
runs along the eastern edge of the galaxy, opposite the 
extraplanar H I gas. Furthermore, the radio continuum 
emission itself appears asymmetric; the radio disk has 
a sharp cutoff along the eastern edge and extraplanar 
emission past the western edge similar to that of the H I 
gas. This is consistent with the picture in which the ICM 
is acting on the eastern (leading) edge of the galaxy disk, 
where we find the radio continuum to be deficient, and is 
pushing the ISM (gaseous and relativistic phases) toward 
and beyond the western edge. 

5.2. Radio Polarization and Spectral Index Maps 

Radio polarization and spectral index data of 
NGC 4522 provide evidence for an ongoing ICM-ISM 
interaction (|Vollmer et al. ir2004a[) : the eastern edge of 
NGC 4522 is found to be highly polarized and has the 
flattest spectral index, which steepens as one moves to- 
wards the western side of the galaxy. In the first column 
of Figure [5] we plot the radio continuum deficit regions 
and overlay the 20 to 6 cm spec tral index and 6 cm po - 
larized radio continuum maps of lVollmer et al. 1 ([2004a) . 
We find that the regions of high polarization and flat- 
test spectral index (i.e. a ~ 0.7) lie just interior to 
the radio deficit region. While the peak in the polar- 
ized radio continuum is coincident with the flattest part 
of the spectral index distribution, this is not where the 
total radio intensity ( or FIR) peaks as pointed out by 
IVollmer et al. I (|2004al ). This is also evident in Figure |H] 
where we show the 6 cm radio continuum emission to- 
gether with the vector s of the magnetic field uncorrected 
for Faraday rotation (jVollmer et al. Il2~004a| ). on the ra- 
dio deficit map. The magnetic field vectors appear to 
be elongated parallel to the radio deficit region and the 
ICM-ISM interaction surface. 

These observations suggest that, except for syn- 
chrotron tails or the radio deficit regions, to first or- 
der, star formation in the disk still drives the appear- 
ance of the FIR/radio ratios. However, by inspecting 
the observed-to-modeled radio continuum ratio map of 
NGC 4522 displayed in the second column of Figure 
[5j which was created using our smoothed 70 /im image, 
ICM-ISM effects become more apparent. We find a lo- 
calized peak in the ratios just outside the radio deficit 
region, and near the regions of high radio polarization 
and flat spectral indices; we attribute this peak to a lo- 
cal enhancement in the radio brightness. We now try to 
determine the most plausible physical scenario to explain 
this combination of observations. 

5.3. Ram Pressure and the Relativistic ISM: Magnetic 
Field Compression, Shearing, and/or CR 
Re- acceleration? 

Ram pressure from the ICM can in principle work 
to both com press the ISM, possibly tr i ggering star for- 
matio n (e.g. iBothun fc Dressier"! 119861 : ICaldwell et al. I 
1993), as well as truncate a galaxy's star-forming 
disk by sweeping away interstellar material (e.g. 
iKennev fc Koopmann I Il999f ). Similarly, ram pressure 
may act on a galaxy's relativistic ISM by compressing 
or stretching its magnetic field, redistributing CR parti- 
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Fig. 4. — Ratio maps of the observed to expected radio flux for six galaxies observed as part of SINGS (see i|3.3l l. The input radio and 
infrared images were smoothed to simulate the expected appearance of each galaxy at the distance of Virgo. 



cles, as well as by re-accelerating CR particles via ICM 
triggered shocks. 

While the high level of radio polarization in NGC 4522 
could arise from a compressed magnetic field, such an 
occurrence is not clearly supported by our observations. 
Since the synchrotron emissivity goes roughly as the 
square of the magnetic field strength, we would expect 
a distinct region of enhanced radio continuum emission, 
where the field is being compressed and amplified, to lie 
just interior to the radio deficit region. If there is such 
an enhancement in the total radio continuum just behind 
the radio deficit region we are probably not resolving it, 
but its amplitude must be modest since a strong signa- 
ture in the ratio map is not detected. 

On the other hand, we also observe a large region of 
relatively flat spectral in dices (i.e. a ~ 0.7) interior to 
the de ficit region. While Nikla sTKlein. fc Wielebinskv I 
(1997) report an average radio spectral index of 0.83 ± 
0.02 for a sample of 74 star-forming galaxies, the radio 
spectral indices of SNRs, whose associated shocks are 
thought to be the primary acceleration sites for CRs, are 
found to be slightly flatter, typi cally ranging between 
~0.5-0.7 both in the Galaxy (e .g. lRevnolds et al. 1120081: 
Allen. Houck. fc Stumer"! I2008D and M 33 (jDuric et al. I 
19951) . These values are also consistent with those 



predicted by diffusive shock acc eleration theory (e.g. 
iDrurv. Markiewicz. fc Volk II 1989( 1. Thus, it seems plau- 
sible that the flat indices interior to the deficit region of 
NGC 4522 is the signature of CR electron acceleration 
sites within NGC 4522. 

We would therefore expect such areas to exhibit mod- 
erate enhancements in total continuum emission. These 
observations are consistent with the interpretation that 
the ICM wind is sweeping up CR electrons in a low pres- 
sure outer-region of the disk as well as driving a series 
of shocks into the ISM which will re-accelerate CR elec- 
trons there. The latter secondary shocks, or shocklets, 



will likely originate at the interface between the ICM and 
ISM, arising from inhomogeneities and turbulence. This 
would then explain the large region of flat spectral radio 
indices and area of modest increase to the radio contin- 
uum, as observed. Some of the CR electrons, including 
those that are newly accelerated, will be swept behind 
the galaxy creating a synchrotron tail that limit the en- 
hancement to the radio surface brightness even though 
acceleration is occurring. 

A second explanation for the existence of the high ra- 
dio polarization could be that of shear as discussed by 
Vollm er et al. I (|2004al ) which does not require total mag- 
netic field amplification. In this scenario the ICM wind 
sweeps up CRs and the magnetic field; the ISM magnetic 
field lines are stretched thereby increasing the fraction of 
polarized radio continuum emission without significantly 
increasing the total field strength. This process alone, 
however, does not explain the observed flat spectral in- 
dices. 

We therefore prefer the scenario in which the magnetic 
field of NGC 4522 is being sheared along the interacting 
region between its relativistic disk and the ICM as ram 
pressure induced shocks propagate into its ISM. The lo- 
cally depressed radio continuum emission in the deficit 
region is the result of ram pressure from the ICM sweep- 
ing up low density CR electrons, associated with past star 
formation activity which have diffused significantly far 
away (~kpc) from their acceleration sites in supernovae. 
While we do not find significant enhancements to the to- 
tal radio continuum behind the deficit regions, the mag- 
netic field may still be moderately compressed. Some of 
the CR electrons which have been swept-up by the ICM 
are redistributed downstream producing the synchrotron 
emission tails as ob served very promine ntly here, as well 
as for NGC 4402 (jCrowl et al. I [2005h ; this effect will 
limit the enhancement to the radio surface brightness. 
Another subset of particles, however, are re-accelerated 
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Fig. 5. — In the top panels, from left to right, we plot the 
radio continuum deficit map and the ratio map of the observed 
to modeled radio continuum emission of NGC 4522 overlaid wit h 
radio spectral index contours taken from [Vollmer et al. I {2004a ). 
The contour levels increase by 0.15 and range between 0.7 to 1.9; 
the dashes indicate the direction of the downward gradient. We 
plot the same images in the bottom panels exce pt this time we 
overla y them with 6 cm polarized radio contours (Vollm er et al. I 
l2004aT l, The contours are logarithmically scaled and begin at the 
3-cr RMS level. The cross in each panel identifies the center of the 
galaxy. 
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Fig. 



- The radio deficit map of NGC 4522 plotted with the 
6 cm total continuum emissions as contours and th e vectors of the 
magne tic field uncorrected for Faraday rotation (| Vollmer et al. I 
2004aj). The contours arc logarithmically scaled and begin at the 3- 
a RMS level. The size of the magnetic field vectors are proportional 
to the intensity of the polarized radio emission. The beam (19" X 
19") is plotted in the bottom left corner of the image. 

interior to the radio continuum deficit region raising the 
total radio power as evidenced by the flat spectral indices 
and modest increase in the radio continuum found in the 
ratio map of NGC 4522. The preference for this scenario 
is discussed further in the context of the observations 
from the entire sample in §6.2. 

6. DISCUSSION 

We have shown that a large fraction of our sample 
of cluster galaxies exhibits statistically significant ra- 
dio deficit regions relative to the prediction of the phe- 



nomenological image-smearing model of iMurphv et al71 
( 2008). Such deficit regions are not found within normal 
field galaxies, suggesting that cluster environment pro- 
cesses are at work. Comparing our results with other ob- 
servational data for NGC 4522, a galaxy which is clearly 
experiencing the effects of ram pressure, we find our 
deficit region agrees with such a scenario and adds new 
insight into the interaction. 

Similar to the example of NGC 4522, radio polarization 
observations of other Virgo cluster galaxies exhibit highly 
polarized radio emission along the side of the galaxy pre- 
sumed to be the leading edge of the ICM-ISM intera ction. 
(|Vollmer et al. I [20071 : IChvzv. Ehle. & Beckl [2007m . For 
the galaxies in our sample having polarization data and 
which show evidence for being disturbed (i.e. NGC 4254, 
NGC 4388 and NGC 4402), the regions of highly po- 
larized continuum emission appear coincident with the 
modest local radio continuum enhancements found in our 
ratio maps (i.e. ratios of ~1.1 — 1.7) and interior to our 
radio deficit regions. When averaged over these kpc-sized 
regions, such enhancements in the total radio continuum 
translate, roughly, into a compression and amplification 
of the magnetic field by ~5 — 30%, clearly a modest ef- 
fect. Thus, the scenario described for NGC 4522 seems 
to be applicable for these galaxies as well based on the 
polarized radio continuum results. With this picture we 
will use the findings from our analysis to quantitatively 
assess the strength of the ram pressure from the intr- 
acluster medium (ICM). Since the radio deficit regions 
detected for NGC 4396 and NGC 4579 differ empirically 
from what we find for the rest of the sample, and show 
evidence for being artificial, we exclude them from the 
following discussion (see ij4.ip . 

6.1. Estimates of Minimum Ram Pressure and 
Internal Relativistic ISM Pressure 

Ram pressure is expected to vary by several or- 
ders of magnitude as galaxies orbit within clusters 
(Voll mer et al. I [20011 : iTonnesen. Bryan, fc van Gorkoml 
l2007f) . but there is no known way to measure the strength 
of current ram pressure on a galaxy. The radio deficit of- 
fers the potential to do this. 

Figure [7J shows that the strength of the radio deficit is 
closely related to the time since peak pressure as inferred 
from st ellar population studies and gas stripping simu- 
lations dCrowl fc Kennev 1120081 IVollmer fc Huchtmeier I 
120031: IVollmer et al. Il2004bl 120061 ). This strongly sug- 
gests that the strength of the radio deficit is good in- 
dicator of the stren gth of the current ra m pressure. In 
simple models (e.g. IVollmer et al. I l200lh the ram pres- 
sure drops by roughly an order of magnitude 300 Myr 
after peak pressure, and by a factor of ~5 from 100 Myr 
to 300 Myr after peak pressure. As shown in Figure [7J 
this is similar to the change in the radio deficit parameter 
between NGC 4522 and NGC 4569, suggesting that dif- 
ferences from galaxy to galaxy in the radio deficit param- 
eter are reflecting differences in the strength of current 
ram pressure. 

We can get a rough estimate of the strength of ram 
pressure acting on these galaxies from the equipartion 
energy density of the radio-emitting plasma. Ram pres- 
sure from the ICM is simply defined as, 
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TABLE 5 

Radio Continuum Deficit Regions Results 



cobs cmod A pmod pobsdisk j-a 

°v °v ^def ^RISM ^RISM 'quench 

Galaxy (mjy) (mjy) X (kpc 2 ) A dct /A disk (HT 12 dyn cm" 2 ) (10~ 12 dyn cm~ 2 ) P» s d u /Pg sk (Myr) 



Virgo (Cluster) Sample 



n4254 


1.578 


4.728 


0.007 


33.0 


0.04 


3.9 


8.3 


0.46 




n4330 


1.891 


4.939 


0.153 


20.7 


0.27 


2.0 


2.1 


0.97 




n4388 


1.733 


7.607 


0.036 


22.6 


0.21 


2.1 


5.6 


0.38 


225 


n4402 


1.120 


4.134 


0.045 


18.8 


0.13 


2.3 


4.3 


0.53 


200 


n4522 


1.676 


4.228 


0.096 


18.3 


0.17 


2.5 


3.2 


0.80 


100 


n4569 


1.852 


4.702 


0.023 


35.6 


0.12 


2.6 


7.1 


0.36 


300 



Field (Control) Sample 

NGC 3627 2.207 7.318 0.011 43.3 0.08 2.9 6.8 0.43 

Note. — NGC 4321 and NGC 4580, along with nearly every field galaxy, are excluded since radio continuum deficit regions were not 
found within their disks; we note that the ^RisM Sk va l ues f° r these two galaxies are 7. 1 an d 16.0 x 10~ 12 dyn cm" 2 , respectively. NGC 4396 
and NGC 4579 were excluded since their radio deficit regions ma y be artificial ( sec § |4. H . 

a : Time since outer-disk star formation was quenched taken from lCrowl fc Kennevl J2008D ; for NGC 4580 they report t qucnch = 475 Myr. 



where picm is the ICM mass density and w ga i is velocity of 
a galaxy relative to the ICM. If the ICM ram pressure ex- 
ceeds the internal pressure of a galaxy's relativistic ISM 
(CRs + magnetic held), then it should be possible to 
redistribute and even strip them from the galaxy disk. 
Using the predicted radio flux density for each deficit re- 
gion, we can approximate a minimum value for the ICM 
ram pressure needed to cause the observed depression in 
the radio. 

Taking the predicted flux density of the deficit region 
along with its area we use th e revised equipar t ition and 
minimum energy formulas of lBeck fc Krause I (|2005f ) to 
calculate the minimum energy magnetic field strength 
of the deficit regions before they were stripped. This 
calculation assumes a proton-to-electron number density 
ratio, for particles in the energy range corresponding to 
GHz synchrotron emission, of 100, a radio spectral index 
of 0.8, and a path length through the emitting medium 
of X/cos(i) kpc where i is the galaxy incli nation. The 
inclina tion is derived using the method of iDale et al. I 
(|1997f) such that, 



cos 2 i 



{b/af - (6/a)? nt 



) 2 

Ant 



(7) 



1 - (b/a) 

where a and b are the observed semi-major and semi- 
minor axes and the disks are oblate spheroids with 
an intrinsic axial ratio (6/ a )int — 0.2 for morphologi- 
cal types earlier than Sbc and (fo/a)i n t ~ 0.13 other- 
wise. Using these calculated magneti c field strengths 
B (i.e . calculated using Equation 4 of iBeck fc Krausel 
(2005)) we compute the magnetic field energy densities 
Ub = B 2 /(8tt) of the deficit regions. 

Assuming minimum energy between the magnetic field 
and CR particle energy densities, Ub and Uqr respec- 
tively, we can use the values of Ub over the deficit re- 
gion in our model to determine the minimum P ram nec- 
essary to create the deficit regions. From minimum 
energy arguments we find that Ub = 3/4<7cr- For 
a relativistic gas the magnetic pressure, Pb, and CR 
pressure Pqb, are related to energy density such that 
Pb = Ub and Pgr = 1/3(7cr leading to the relation 
that Pgr = 4/9E/b- Then, the pressure of the relativistic 
ISM is found to be 

Prism = Pgr + Pb - 13/W B . (8) 
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Fig 7. — The severity of the deficit region characterized by T 
(see i|3,3H plotted a gainst time since star fo rmation was quenched 
in their outer disks (Crowl & Kcnney 2008), likely coinciding with 
the time since the galaxy disk experienced peak-pressure. 
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Fig 8. — The severity of the deficit region characterized by T 
(see i|3.3l l plotted against a measure of TIR surface brightness. 



The relativistic ISM pressure of the deficit regions before 
they were stripped, as estimated from our models and 
denoted as P£g M , are given in Table these estimates 
set the minimum ICM ram pressure necessary to create 
the observed radio continuum deficit regions (i.e. P ra m > 

pmod \ 

All of the Virgo galaxies listed in Table [5] have Prism 
values ranging between ~2 — 4 x 10~ 12 dyn cm~ 2 . Us- 
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ing the projected linear distances to the cluster center for 
these sample galaxies along with a measured density pro- 
file of Virgo ([Matsumoto et al. 1 1200(1 yields a range m 
ICM density of nicM = Picm/"i p « 1.3-4.2xl0~ 4 cm~ 3 . 
Taking a Virgo galaxy velocity of 1000 km s _1 , which 
is comparable to the ~930 km s" 1 velocity dispersion 
amon g Virgo spirals ( Bin ggeli. Tammann. fc Sandage I 
1 1987T) and typically achieved by cluster members dur - 
ing some portion of their orbits (|Vollmer et al. I 2001). 
the associated range in ICM ram pressure is P ram « 
2 — 7 x 10~ 12 dyn cm -2 . This range of ICM pressures is 
generally greater than or similar to the values for P^j^ 
which is consistent with ICM ram pressure being able 
to create the observed deficit regions. These pressures 
agree with t hose produced in the 3D h ydrodynamical 
simulation of iRoediger fe Briiggen I (|2007f ) over the same 
projected linear distances to the cluster center among 
our sample. We also note that, given the ICM density 
range, the corresponding velocities required to reach a 
minimum ram pressure of ~2 x 10~ 12 dyn cm~ 2 range 
between 530 and 980 km s , quite comparable to the 
velocity dispersion of Virgo spirals. 

While this value only sets the minimum ram pressure 
strength required to create the deficit region, it is per- 
haps more interesting to try to get a better handle on 
the equilibrium pressure at the working surface between 
the deficit region and the edge of the galaxy's observed 
radio disk. Since it is difficult, and probably not very 
accurate, to calculate the relativistic ISM pressure on 
a point-by-point basis within the disk, we instead mea- 
sure the average internal relativistic plasma pressure of 
observed radio disk P^|^ sk . We define P^|^ sk using 
the observed total radio flux density and the area given 
within the 3-cr RMS isophote of each radio image. The 
values of -PpigM Sk : along with the ratios between Pr^sm 
and PjS s \ are given in Tabled While NGC 4321 and 
NGC 4580 were excluded in Table [5] since they do not 
have radio deficit regions, we note that their values of 
^R.iSM Sk are 7-1 an d 16.0 x 10~ 12 dyn cm -2 , respectively. 
Therefore, the full range of relativistic plasma pressures 
are similar to ram pressures calculated above for a typical 
Virgo galaxy and a range of ICM densities. 

The internal relativistic plasma pressures cannot be 
vastly different from that of the ICM ram pressure since 
then we would either not observe the radio deficit regions 
or the entire relativistic ISM would be swept out of the 
disk. We expect to see more significant radio deficit re- 
gions (i.e. larger values of T) for cases where the internal 
relativistic ISM pressure is low. Comparing the values of 
T with -P£^ sk in Table we find that this is indeed 
the case. Our ram pressure estimates from minimum en- 
ergy arguments give a good rough approximation for the 
strength of ram pressure within a factor of a few, but can- 
not be expected to measure precisely the differences in 
ram pressure among the sample galaxies. Further work 
is required to more closely and quantitatively relate the 
radio deficit to the strength of ram pressure. 

In Figure [8] we find that there is a general trend of de- 
creasing values of T with increasing TIR surface bright- 
nesses as measured by the ratio of each galaxy's global 
TIR flux to Adisk, the deprojected area of their ratio map. 
Since TIR surface brightness and pRigM Sk are directly re- 
lated, this trend is not terribly surprising. Interpreting 



the TIR surface brightness to be a proxy for the mean 
ISM density, we find that galaxies having a more dense 
ISM appear to be less affected by the ICM. This seems 
reasonable given that the ram pressure would need to 
be stronger to disturb their ISM and is consistent with 
higher values for the internal relativistic plasma pressure 
for small values of T. However, this trend may be for- 
tuitous given that the external ram pressure likely varies 
tremendously among the sample. 

6.2. Combining Local and Global Observations into a 
Single Physical Scenario 

Using our results on the spatially resolved FIR/radio 
distributions of these cluster spirals, along with their 
global FIR/radio properties, we will now investigate 
whether the ICM wind plays a significant role in enhanc- 
ing the escape and re-acceleration of CR electrons. Af- 
ter considering escape, we explore a number of physical 
scenarios to explain simultaneously the local and global 
FIR/radio ratio behavior among these galaxies before 
coming to our preferred picture. 

6.2.1. Cosmic Ray Electron Escape 

The only galaxy with a high global q value is 
NGC 4580. While the CR electrons are being moved 
around by the ICM wind, it is unclear as to whether 
they easily escape the galaxy disks as global q (logarith- 
mic FIR/radio) ratios do not appear systematically high 
wit h respect to the nominal value of ~2.34 ± 0.26 dex 
(i.e. lYun. Reddv fe Condon 112001f t. In fact, by plotting 
q against T (our parameter defining the severity of the 
deficit region) in Figure O we find a trend of decreasing 
FIR/radio ratios with increasing T and that nearly all 
of the q values are lower than the nominal value. Other, 
more detailed studies on the global FIR and radio prop- 
erties of cluster galaxies have pointed out notably lower 
FIR/radio ratios t han expected from the FIR-radio re- 
lation in the field ([Miller fc Owen 1 120011 : iReddv fc Yu"n1 
2004). This suggests that there is either depressed FIR 
emission or an enhancement in the global radio emission 
among our sample galaxies. And while the trend in Fig- 
ure[5]may simply be the result of small number statistics, 
it also might suggest that q is sensitive to the strength of 
the ICM ram pressure and perhaps a galaxy's stripping 
history. 

We do not detect a radio deficit region for NGC 4580. 
This galaxy is highly H 1 deficient, having the largest 
deficiency value in Table Q] (1-53 dex; a factor of ^5 
larger than the median value among t he sample), and 
has b een classified as an anemic spiral ([van den Bergh I 
119761) . though better de scribed as having a severely tr un- 
cated star- forming disk (Koopm ann fc Kennev Il2004h . It 
appears to have been stripped long ago (^475 Myr) and 
is the only post-strongly stripped gal axy in our sample 
obser ved so long after peak pressure ([Crowl fc Kenneyl 
2008) . We therefore speculate that its high q value is the 
result of CR electron escape as the ICM wind has swept 
out most of the galaxy's gaseous and relativistic ISM 
long ago leaving only the central disk. The only other 
galaxy for which we do not detect a radio deficit region 
is NGC 4321 which is much different from NGC 4580. 
NGC 4321 has not been strongly stripped and has a near- 
average q value; this galaxy appears much more similar 
to the WSRT-SINGS field galaxies. 
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6.2.2. Low Global FIR/Radio ratios as a Result of 
Re-acceleration? 

As already discussed in £15.31 we have attributed local 
enhancements to the radio continuum arising from the 
re-acceleration of CR electrons by ICM induced shock- 
lets. While moderate local enhancements are found 
in the radio continuum interior to the radio deficit re- 
gions via the ratio maps, they do not appear signif- 
icant enough to explain q ratios which are <0.4 dex 
lower than average (i.e. a factor of <2.5). The q val- 
ues of 'Taffy' systems are found to be a factor of ^2 
low due to excess radio continuum emission associated 
with the synchrotron bridges connecting the galaxy pair s 
(jCondon et al I Il993t ICondon. Helou. fc Jarrettl l2002f ). 
Looking to see if synchrotron tails may be causing a sim- 
ilar effect in NGC 4522 (where the observed synchrotron 
tail appears most striking) we find that this cannot be 
the case. Only 15 — 20% of the total radio continuum 
emission is associated with the tail. 

While we have argued that it does not seem likely that 
the magnetic field has been compressed significantly (i.e. 
no more than <30% as suggested by the radio excesses) 
at the leading edge between the ISM and the ICM wind, 
leaving no clear signatures of a gradient in the ratio 
map, perhaps the field strength is raised globally due 
to therma l compression by the su rrounding ICM as sug- 
gested bv lMiller fe Owen I (|2001f ). Models of such an ef- 
fect seem to be able to produce the factor of 2 — 3 excess 
sync hrotron emission and reproduce the observed q ra- 
tios (|Reddv fc Yun |[2~004h . If such a scenario were true, 
and the Virgo galaxies with low q ratios are simply com- 
pressed versions of field galaxies, we would expect galax- 
ies with the lowest q ratios to appear compressed in all 
directions and exhibit relatively larger internal relativis- 
tic plasma pressures, ^RisM Sk ' compared to Virgo galax- 
ies with larger q ratios. However, the galaxies having the 
lowest q ratios both exhibit extraplanar radio continuum 
emission along their trailing edges in the form of diffuse 
synchrotron tails and do not have larger -Prism 515 values 
relative to the rest of sample, inconsistent if compression 
was the dominant mechanism. 

Another possibility is that the FIR emission is de- 
pressed in these systems through the removal of a dust 
component. For instance, the ICM wind may have re- 
moved a substantial amount of the cool dust (cirrus) 
component, a significant contributor to the total FIR 
emission. However, FIR color differences between the 
stripped Virgo galaxies and the WSRT-SINGS spirals 
are too small to induce differences in the SEDs and to- 
tal FIR emission large enough to account for the dif- 
ferences in q. The MIPS /„(70 //m)//„(160 fim) and 
/„(24 yum)// !/ (70 /xm) flux density ratios have medians 
that are < 20% lower and higher, respectively, than those 
found for the WSRT-SINGS spirals. Both samples have 
FIR colors in the range of quiescent normal star-forming 
galaxies, where the 20% color differences translate into 
small changes in the SED, and changes of only a few 
percent in FIR norma lized either to t he 24 or 70 /im 
emission, based on the iDale fc Heloul (|2002l ) SED tem- 
plates. The 24 /im emission component is relevant for 
normalization because of its close assoc iation with HII 
regio ns and dense, star-forming regions (Cal zetti et al. I 
l2007f ). making it much more resistant to stripping than 



the longer wavelength c omponents. Th i s resu lt is con- 
sistent with the study of iReddv fc Yun I (|2004h who also 
found that the FIR properties (i.e. luminosities and col- 
ors) of cluster galaxies having excess radio emission are 
similar to those of field galaxies. 

The FIR emission in these systems may also be de- 
pressed, without systematic changes in their FIR col- 
ors, as a result of depressed star formation activity per 
unit area. If ram pressure from the ICM wind is driv- 
ing shocklets into the ISM, mechanical energy would be 
dissipated into molecular clouds as turbulence, possibly 
reducing their efficiency to collapse and form new stars. 
However, in order for a galaxy to exhibit q values simi- 
lar to those reported here, the decrease in star formation 
activity must have begun within the last ^50 Myr (i.e. 
roughly the expected CR electron cooling time to syn- 
chrotron losses). At later epochs, the synchrotron emis- 
sion associated with the previous star formation episodes 
fades away, and q values drift back towards more nor- 
mal values even if the star formation activity remains 
depressed. Since 50 Myr is more than an order of mag- 
nitude smaller than the cluster crossing time, it is un- 
likely that we would be observing so many galaxies dur- 
ing this special period of time. It therefore appears that 
depressed FIR emission is not the dominant effect driving 
the systematically low q values for these galaxies. 

While it is still not completely unambigous, we be- 
lieve the most plausible explanation is that the ICM 
wind raises the global CR electron energy and syn- 
chrotron power by driving shocklets into the ISM which 
are re-accelerating CR electrons. Taking again i> ga i sw 
1000 km s" 1 and the range in picm values discussed in § 
16. D and assuming the average surface area of each galaxy 
disk exposed to the ICM wind is 25 kpc 2 (i.e. the mean 
size of the radio deficit regions given in Table [5]), we com- 
pute the kinetic energy per unit time a galaxy will inherit 
from the ICM wind. Comparing this number to that of 
the mean radio continuum luminosity among those galax- 
ies listed in Table[5j and assuming a number density ratio 
of protons-to-electrons, for particles in the energy range 
corresponding to GHz synchrotron emission, of 100, we 
find that ^3 — 11% of the total kinetic energy available 
must be given to the CRs in order to double the radio 
luminosity While a value ~11% may be slightly large, 
there are many uncertainties involved with this calcula- 
tion including whether the CR electrons are more effi- 
cient than CR nuclei at absorbing this energy. In any 
case, it is certainly clear that the available mechanical 
luminosity of the system is enough to cause the observed 
global excess in synchrotron emission. 

The fact that we only find relatively moderate local 
enhancements to the total continuum behind the radio 
deficit region, save perhaps NGC 4254 where a less mod- 
est radio excess is observed (i.e. ratios peaking at ~1.7), 
suggests that the shocks run through the entire galaxy 
disk rather quickly. The shock speed through the galaxy, 
v s , which must be super- Alfvenic to re- accelerate CR 
electrons, should have a value around 

v s ~ r « Wgal , (9) 

V PISM / V O piSM J 

where P ram is the ram pressure, as defined in Equation 
O and pism and picm are the ISM and ICM densities, re- 
spectively. Again taking v ga \ m 1000 km s _1 , the range in 
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Fig. 9. — The severity of the deficit region characterized by T 
(see £|3.3| I plotted against q (the logarithmic FIR/radio ratio). The 
vertical line at q = 2,34 i dentifi es the average q value reported by 
lYun, Reddv, fc Condonl (2001) for 1809 galaxies. 
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Fig. 10. — The severity of the deficit region characterized by 
T (see £|3.3| I plotted against the distance, in degrees, to the Virgo 
cluster center. 

Picm values discussed in § 16.11 and nisM = Pism/«i p ~ 
1 cm -3 leads to shock velocities ranging between ~15 
to 25 km/s. Assuming a thin disk thickness of 500 pc, 
the shocks should run through each disk within ~20 to 
40 Myr; indeed this is very short compared to the dy- 
namical timescale of these systems. 

As previously noted, this proposed mechanism may 
work to slow down star formation itself, decreasing the 
infrared surface brightness, either from an increase in 
cold molecular cloud heating by the re-accelerated CR 
particles or similar gas-phase shocklets perturbing the 
clouds. However, it is unclear how significant this effect 
may be relative to the increase in synchrotron emissivity 
via the associated re-acceleration of CR electrons. 

6.3. Pre- or Past- Peak Pressure? 

Since many cluster galaxies have highly radial orbits, 
the direction to the cluster center can help provide some 
insight as to whether a galaxy may be traveling into 
or away from the cluster center. The direction toward 
the cluster center is indicated on the radio deficit maps 
in Figure [TJ Some galaxies have deficit regions facing 
toward the cluster center (NGC 4254, and NGC 4402) 
which is suggestive of a galaxy entering the cluster and 
at a phase of pre-peak pressure. On the other hand, we 
also find an example for which the radio deficit region 
faces away from the cluster center (NGC 4569) which 
may suggest that the galaxy is on their way out of the 



cluster and past peak-pressure. There are also galaxies 
which face neither directly towards or away from the clus- 
ter center (NGC 4330 and NGC 4522) leading to a more 
complex scenario; such objects may be stripped by inter- 
actions with sub-cluster structures (e.g. NGC 4522 with 
M 49). We also note that NGC 4388 may be interacting 
with the M 86 sub-cluster; its H I tail points neither to- 
wards or away from the cluster center but rather towards 
M 86 (jOsterloo fe van Gorkoml[200l . 

Another interesting case is that of NGC 4254 which is 
thought to be entering Virgo for the first time at high 
speed . It has recently been suggested b y iKent et al.~l 
(|2007f) and lHavnes. Giovanelli. fc Kent I (|2007t ) that a 
large H I tail associated with NGC 4254 is a tidal arm 
resulting from galaxy haras sment (|Moore et al. I [19961 
iMoore. Lake, fc Katz I I1998D . However , no interacting 
counterpart has been identified. Also, the H I tail is lo- 
cated opposite the observed radio deficit region and the 
direction to the cluster center. Since it appears that the 
radio deficit regions are related to the effects of ram pres- 
sure, it therefore seems plausible that the large H I tail 
could be the result of ram pressure stripping. However, 
it is also possible that the radio deficit region may be 
created via galaxy harassment or that, more plausibly, 
both physical processes are currently at work. 

Projection effects, specifically the orientation of the 
galaxy disk with respect to the ICM wind direction and 
our line-of-sight, likely play a significant role in determin- 
ing the location of deficit region centroids. For largely 
face-on systems, the relationship between the deficit re- 
gion location and the direction of the wind is probably 
more easily predicted than for more highly inclined galax- 
ies. The highly inclined galaxies NGC 4330, NGC 4402, 
and NGC 4522, have radio deficit regions consistent with 
lunes along one side of their disks; these could occur over 
a large range of wind angles. Furthermore, the size of the 
observed deficit region is subject to projection effects as 
well, likely appearing much smaller for disks viewed face- 
on, especially if they are moving through the ICM along 
the line-of-sight. 

Therefore it is hard in general to gauge whether galax- 
ies may be pre- or past-peak pressure from the direction 
of the cluster center with respect to the deficit region. In 
Figure fTUl we plot T against the distance to the cluster 
center. The lack of observed trend suggests that the 3D 
structure of the cluster along with presence of sub-cluster 
structures are important in understanding ICM-ISM in- 
teractions and the picture of a smoothly distributed and 
static ICM is overly simplistic for characterizing ICM- 
ISM interactions. 

6.4. Departures from Typical ISM Conditions 

While we can explain regions having radio continuum 
emission without corresponding FIR emission as the re- 
sult of CR electrons being swept into a synchrotron tail 
by the ICM, the observation of FIR emission without 
corresponding radio emission is a bit more interesting 
since it suggests a strong departure from typical ISM 
conditions. Even though such regions are close to the 
noise limits of the infrared and radio maps, this result 
appears to be significant. In Figure [5] we show two ex- 
amples (i.e. NGC 4402 and NGC 4522) for which we 
detect 70 /jm emission more radially extended than the 
radio; these regions are found only along the galaxy edge 
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associated with the radio deficit region. Furthermore, if 
q were constant across these galaxy disks, we would ex- 
pect to detect the radio continuum in these regions at 
>3-<7 RMS level of the radio maps. However, we know 
that q does not remain constant across galaxy disks and 
actually decreases wit h increasing galactocentric radius 
iMurphv et al. I (|2006aD , implying even higher levels of ra- 
dio continuum should be expected, easily exceeding >6-cr 
RMS noise levels. We therefore conclude that the non- 
detection of radio emission is significant at least at the 
~6-cr level or, in other words, <15% of the expected radio 
continuum emission is present. 

These observations then suggest that the diffuse rela- 
tivistic and gaseous ISM have been stripped away along 
these leading edges leaving the denser molecular mate- 
rial as traced by the 70 /jm dust emission. Presumably, 
the low density relativistic phase is quickly removed due 
to an increase in open magnetic field lines leading to the 
quick escape of CR particles. While our observations do 
not rule out that synchrotron emission associated with 
the remaining molecular clouds is still present and sim- 
ply below our detection limit, since the diffuse compo- 
nent will dominate the total radio continuum emission 
in these regions, it is still interesting to speculate how 
quickly this synchrotron emission associated with the re- 
maining molecular clouds would fade as CR electrons 
diffuse away. 

Assuming a flux freezing scaling of B oc n 1 / 2 (e.g 



Ruzmaikin. Sokoloff fe Shukurov lfl988t iNiklas fe Beck I 
19971: ICrutcher II1999D . a solar neighborhood ISM density 
and magnetic fiel d strength of rtRN = Psn / rn v ~ 1 cm" 3 
and -Bsn ~6(iG (|Strong. Moskalenko. fc Reimer IfeOOOl) , 
and an average giant molecular cloud (GMC) density of 
n-GMC = PGMc/mp w 10 3 cm 
netic field strength of -Bgmc 
actually quite consistent with field strengths of molecu- 
lar clo uds inferred from the Chandrase khar-Fermi tech- 
nique (|Chandrasekhar fc Fermi I I1953T) using polarize d 
dust continuum emission (e g ICrutcher et al. I I2004D . 
Rewriting Equation 24 from IMurphv et al. I ( 2008h for 
diffusion alone, such that 



3 leads to a GMC mag- 
190 /xG. This value is 



Tdiff * = 2 x 10 6 
yr 



kpc 



( " 



GHz 



-1/4 



B 



1/4 



. ( 10 ) 

where t^is is the time it takes a CR electron emitting 
at frequency v to diffuse a distance Zdiff in a magnetic 
field of strength B assuming random walk diffusion char- 
acterized by an energy-dependent diffusion coefficient, 
we can calculate the time it would take the CR elec- 
trons to diffuse away from the molecular clouds. Taking 
B = Bgmc ~ 190 ^G and idiff ~ 0.05 kpc leads to a dif- 
fusion time of ^1.7xl0 4 yr for CR electrons to leave the 
molecular clouds; this is considerably short compared to 
the stripping timescales and plausible given our observa- 
tions. 

Such a scenario may have implications for star forma- 
tion in these regions as well. According to Equation [8] 
and assuming minimum energy arguments, the complete 
lack of CR pressure, which works to support a cloud from 
collapse, will decrease the total relativistic ISM pressure 
by ~30%. This situation may encourage cloud collapse 
in these regions leading to enhanced star formation. On 
the other hand, the ICM wind will also perturb molecular 



clouds introducing additional turbulence further discour- 
aging collapse. Determining which of these competing 
effects will actually dominate is beyond the scope of this 
paper. 

7. CONCLUSIONS 

We have studied the interstellar medium (ISM) of 10 
Virgo galaxies included in the VLA Imaging of Virgo 
in Atomic Gas (VIVA) survey using Spitzer MIPS and 
VLA 20 cm imagery. By comparing the observed ra- 
dio continuum images with modeled distributions (i.e. 
appropriately smoothed FIR images), created using a 
phenomenologi c al image-sme aring model described by 
IMurphv et aLl (|2006al 120081 ), we find that the edges of 
many cluster galaxy disks are significantly radio defi- 
cient. These radio deficit regions are consistent with be- 
ing leading-edge regions affected by intracluster medium 
(ICM) induced ram pressure as suggested by the location 
of H I and radio continuum tails. From our results we 
are able to conclude the following: 

1. The distributions of radio/FIR ratios within cluster 
galaxies thought to be experiencing ICM-ISM ef- 
fects are systematically different from the distribu- 
tions in field galaxies. Radio/FIR ratios are found 
to be systematically low, due to a deficit of radio 
emission, along the side of the galaxy experiencing 
the ICM wind. 

2. In the case of NGC 4522, a clear example of ongo- 
ing strong stripping, we find that the radio deficit 
region lies exterior to a region of high radio polar- 
ization and a flat radio spectral index. We inter- 
pret this to suggest that CR electrons in the halos 
of galaxies are being swept up by the ICM wind. 
The ICM wind drives shocklets into the ISM of 
the galaxy which re-accelerate CR particles inte- 
rior to the working surface at the ICM-ISM inter- 
face. Some CR electrons may also be redistributed 
downstream creating synchrotron tails as observed 
for NGC 4522. The high radio polarization is prob- 
ably mostly the result of shear as the ICM wind 
stretches the magnetic field; compression may also 
play a role, though a modest one, since the total 
radio continuum in these regions does not appear 
significantly enhanced. 

3. The severity of the current ICM-ISM interaction, 
as measured by the normalized difference between 
the modeled and observed radio flux densities in the 
radio deficit regions, appears to be strongly corre- 
lated with the time since peak pressure as mea- 
sured by stellar population studies and gas strip- 
ping simulations. Thus the radio deficit appears to 
be a good tracer for the strength of current ram 
pressure. 

4. Using the identified radio deficit regions we are 
able to get a quantitative estimate of the mini- 
mum strength of the ICM pressure required to af- 
fect a galaxy disk; we find values in the range of 
^2 — 4 x 10~ 12 dyn cm~ 2 . These pressures are 
generally smaller than, but similar to, those esti- 
mated for typical values of the ICM gas density 
and galaxy velocities, as well as the range of ram 
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pressures calculated by 3D hydro-dynamical simula- 
tions. Therefore, our estimates are consistent with 
the scenario of ram pressure creating the observed 
deficit regions. 

5. The internal relativistic plasma pressures range be- 
tween ~2 and 16 x 10 -12 dyn cm -2 and appear to 
be lowest in some of the galaxies with the strongest 
radio deficits. Thus the strength of the radio deficit 
may also be affected by the internal relativistic 
plasma pressures of the galaxies as well as exter- 
nal ram pressure. Further work will be required to 
learn how much the radio deficits depend on these 
2 variables. 

6. The global radio/FIR ratios of these cluster galax- 
ies are systematically higher than the average value 
found for field galaxies and appear to increase with 
increasing severity of the local radio deficit, which 
may trace current ram pressure. We attribute this 
to an increase in the CR energy density which is 
greater when the interaction between the ICM-ISM 
is more energetic. We also show that the amount 
of mechanical luminosity available from a galaxy's 
motion through the ICM is significantly larger than 
what is required to heat the CR plasma. 



This paper describes a preliminary study comparing 
the FIR and radio emission distributions of Virgo cluster 
galaxies using the first 10 out of ~40 galaxies for which 
Spitzer data was available. We have clearly shown that 
by comparing the FIR and radio continuum distributions 
of these cluster galaxies we can address a number of phys- 
ically important characteristics affecting their evolution 
including the strength and directions of the ICM wind. 
The inclusion of the remaining ~30 sample galaxies along 
with other ancillary data sets will help to better distin- 
guish among scenarios presented here and better illu- 
minate our understanding of how the ICM affects the 
relativistic phase of the ISM. 
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